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THEORY OF GALVANOMAGNETIC EFFECTS IN METALS 

by Gale Fair 

Lewis Research Center 


SUMMARY 

The Boltzmann equation for the distribution function of conduction electrons in a 
metal is used to describe the response of the electronic system to applied magnetic and 
electric fields. The vector mean-free-path formulation, originally developed by Taylor 
to calculate the zero-magnetic-field resistance of metals, is used to find an exact solu- 
tion to the linearized Boltzmann equation for static and homogeneous applied fields. The 
conductivity tensor is evaluated in terms of two operators, one describing the scattering 
between electronic states and the second including the effect of the Lorentz force. The 
eigenfunctions of these operators provide a natural framework for a quantitative descrip- 
tion of the galvanomagnetic properties of a metal, for both an assumed Fermi surface 
and a scattering model for the conduction electrons. 

Of particular importance is the treatment of the electronic scattering processes, a 
treatment sufficiently complete so as to allow a separation of effects due to the scatter- 
ing processes from the effects arising from topological features of the Fermi surface. 

To demonstrate the implementation of the theory, and to illustrate the application of the 
representation theory of finite groups, two examples are considered. The free-electron 
gas, with spherically symmetric scattering, is shown to give the result familiar from a 
relaxation-time approximation (i. e. , a Hall effect but no magnetoresistance). The 
second example is a treatment of a degenerate scattering kernel, still within the frame- 
work of a spherical Fermi surface. A special case of this degenerate kernel theory is 
given in a recent work by Jones and Sondheimer, and a sample calculation is presented 
for comparison with this work. 


THE BOLTZMANN EQUATION 

A galvanomagnetic experiment measures the voltage response to an externally ap- 
plied current for a metal sample in a magnetic field. The experimentally determined 
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quantity is the resistivity tensor of the metal as a function of the magnetic field (both 
strength and orientation) and the temperature. The experiment constitutes a probe of 
the properties of the conduction electrons in the metal and should, in principle, yield 
information about the scattering processes present and the structure of the Fermi sur- 
face. Although galvanomagnetic experiments have been performed for many years, no 
general theory has yet appeared that gives a quantitative connection with experiment. In 
this report a semiclassical theory is presented for the galvanomagnetic effects in metals 
within the framework of the independent-particle description and the Boltzmann equation. 
A closed-form solution of the Boltzmann equation is found that is valid for any strength 
of the magnetic field, and the resulting conductivity tensor is calculated. After a dis- 
cussion of the implementation of the formalism for calculations, consideration is given 
to the special case of a metal having a spherical Fermi surface, for two different scat- 
tering models. 

The semiclassical theory for electrons in metals is based upon the existence of a 
distribution function f(r,k, t). The quantity f(r,k, t)drdk is the probability of finding 
an electron within the element of volume dr dk that is centered at the point (r, k) at the 
time t in the six-dimensional space of position r and Bloch wave-number index k 
(ref. 1). For a spatially homogeneous metallic sample with no thermal gradients pres- 
ent, it is customary to assume that the spatial variation of f is negligible (ref. 2). 

Then for f we can write f(k, t), or f ^ (t) . For a system of electrons in equilibrium, 
the Pauli exclusion principle determines the distribution of the electrons in momentum 
space. The distribution function depends only upon the energy of the electron and is the 

Fermi function, = jl + exp[(e^ - £)/kgT|j"*; where is the energy of the elec- 

tron, T is the temperature of the system, £ is the chemical potential, and kg is the 
Boltzmann constant. This function implies a discontinuity in the energy distribution of 
the electrons when the system is at absolute zero. When T = 0, f*\e k ) is 1 if < £ 
and zero if e k > The value of the chemical potential at absolute zero £(T = 0) de- 
fines the parameter €p, the Fermi energy. For the special case of a system of free 
electrons, the occupied states in wave-number space form a sphere, of radius 
kg = j/ 2mep/h, known as the Fermi sphere. 

The existence of a discontinuity in the distribution function as a function of wave 
vector k is very useful for the discussion of the transport properties of metals. The 
locus of this discontinuity in k- space at T = 0 is called the Fermi surface. Although 
the Fermi surface for the system of conduction electrons in a metal is distorted from a 
sphere by the presence of the lattice of ion cores, experimental evidence (ref. 3) shows 
the existence of such a discontinuity surface at low temperatures. The galvanomagnetic 
experiment is one of a class of experiments designed to probe the Fermi surface of 
metals (refs. 3 and 4). As will be seen later, it is the Fermi surface that helps to de- 
termine the k- states that may participate in transport effects. 
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Before the Boltzmann equation is presented, it is useful to briefly summarize the 
semiclassical view of a conduction electron in a metal. At a finite temperature, the 
discontinuity in the distribution function expressed in terms of is not sharp, the en- 
ergy range which separates the totally occupied states from the totally unoccupied states 
being a band of width approximately kgT. Since this thermal broadening constitutes a 
range of energy which is small compared with the Fermi energy, there is still an ex- 
perimentally well-defined Fermi surface, whose detailed shape is determined by the 
number of valence electrons available, the symmetry of the ionic lattice, and the inter- 
actions between the electrons and their environment (ref. 3). An electron in a state k 
at the Fermi surface moves with the velocity V k e(k)/h, which is, in fact, the group 
velocity of the wave packet describing the electron (ref. 5). An electron in a state k 
may undergo a scattering, which is pictured as a discontinuous change from the initial 
state k to a final state k' . The final state may or may not have the same energy as the 
initial state, depending upon the nature of the scattering process. Common causes of 
scattering (ref. 2) are impurities, physical defects in the crystal such as point defects or 
dislocations, and, at a finite temperature, the thermal vibrations of the lattice (i. e. , 
the phonons). A satisfactory theory of the transport properties of a metal must be 
capable of describing the response of the electronic system to external forces, and in- 
clude specifically the effects due to the scattering processes and the nature of the Fermi 
surface. 

The semiclassical description of transport phenomena is governed by the Boltzmann 
equation (BE) for the electronic distribution function. The BE is the equation of motion 
for the distribution and is obtained by a consideration of the processes which determine 
the net change in occupation of the volume element dr dk. Although the system is not in 
an equilibrium state when external electric and magnetic fields are present, it is cus- 
tomary to assume steady- state conditions. This means that there is no net flow from 
the volume element under observation, and thus the distribution function is not a function 
of time. The effect of the external forces must balance the rate of change of occupation 
caused by the scattering. The equation of continuity, applied to the six-dimensional 
phase space volume element, gives 



= -V £ . *£ 
r k dt 


- V f ^k 
k k * dt 



( 1 ) 


where the subscripts on the gradients indicate the variables of the differentiation and the 
term (3f/8t) co ^^ is symbolic for the rate of change in f caused by scattering. In a gal- 
vanomagnetic experiment the temperature is uniform throughout the sample. In this 
case (ref. 2), the spatial derivative of f k will vanish. For a Bloch electron in electric 
and magnetic fields, the time rate of change of the crystal momentum is given by the 
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Lorontu force law (ref. 6) with the velocity of the electron interpreted as ft” 1 Se/ 9k. 
Thus. 


R ^ = e^E + i v k x (2) 


The scattering term is obtained by considering the balance between scattering into 
ana from the state k. Let Q(k, k’) represent the transition probability per unit time for 
scattering from the initial occupied state k to the final empty state k'. The occupation 
c: the initial and final states is restricted by the Pauli exclusion principle for the elec- 
trons. and so the scattering rate must include factors to ensure that the initial state is 
occupied while the final state is initially unoccupied. The net change in the occupation 
of the state k due to scattering is (ref. 7) 



The volume element in the integral of equation (3) is defined as (ref. 2) 


(2tt) 3 

where V is the volume of the sample. The transition probability per unit time Q(k, k') 
provides the connection with the microscopic description of the electronic properties of 
metals. This probability of scattering from a state k to a state k' is proportional to 
the square of the matrix element of the transition matrix for the scattering between the 
initial and final states (ref. 2). 

The form of the scattering integral (eq. (3)) may be simplified by the assumption 
that the scattering is elastic. Under this assumption the transition probability is sym- 
metric (ref. 2), that is, 

Q(k,k') = Q(k-,k) 

where Q includes a factor 6(e k - to conserve energy between the initial and final 
states. This simplification gives from equation (3), 



f Q(k,k’)[f(k’) - f(k)]dk' 


( 4 ) 
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And the BE for elastic scattering is 


+ - 1 - v k x Hj • V k f(k) = f Q(k, k’)[f(k) - f(k^dk* (5) 

with the energy-conserving 6-function implicit in Q(k, k'). 

There are two questions that have been passed over in the derivation of the BE 
(eq. (5)). The first question is the effect of the omission of inelastic scattering proc- 
esses in passing from equation (3) to (4). The second is concerned with the meaning of 
the integral in equation (4). The symmetry of the transition probability holds true only 
when the scattering is elastic, as assumed. A more general relation may be derived 
from the observation that in the equilibrium state, in which the distribution function is 
the Fermi function, the integral of equation (3) must vanish, so that 

Q(k, k’)f°(e k )[l - f°(e k ,)]= Q(k , ,k)f°(e k ,)[l - f°(e k )] 

Thus, even though Q is not symmetric, there is a function (ref. 7) 

P(k,k’) = Q(k,k’)f°(€ k )[l - f°(e k ,)]= P(k\k) (6) 


that has the desired symmetry. In terms of this function, the time rate of change of f k 
due to arbitrary scattering processes is 



f (k’)[l - f(k)] 

mE 1 - f0 ( e k)] 



( 7 ) 


In principle, the formalism which is presented in this report may be carried through by 
using this more general form for the scattering term. This, in fact, was the approach 
followed by Taylor (refs. 8 and 9) in his solution of the zero-magnetic-field BE. Inelas- 
tic electron-phonon scattering makes a negligible contribution to galvanomagnetic effects 
at the low temperatures at which these experiments are performed (ref. 7). The in- 
clusion of only elastic scattering in the present work should still provide a useful com- 
parison with experiment. 

The problem involved in defining the range of the k- space integral of equation (4) is 
a question of nonuniqueness commonly encountered in the band theory of metals. This 
nonuniqueness has its origins in the translational symmetry of the lattice. In a pure 
metal crystal of infinite extent there are a countably infinite number of points in the 
metal which are equivalent. These points are connected by finite displacement vectors, 
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not by infinitesimal transformations. A translational invariance in configuration space 
always has implications in the momentum properties of a quantum -mechanical system 
(ref. 10). In band theory this gives rise to the definition of the Brillouin zone and the 
use of two pictures, the extended zone scheme and the reduced zone scheme. The ex- 
tended zone scheme permits all values of k; equivalent points on the multiply occurring 
Fermi surfaces are connected by reciprocal lattice vectors. In the reduced zone scheme 
a Brillouin zone is constructed about a selected lattice point in k- space, the Brillouin 
zone being the locus of points nearer this selected lattice point than any other. In this 
case the only scattering transitions allowed are those which connect two points on this 
single Fermi surface (which, of course, may be composed of a number of separate 
pieces). A much clearer description of the reduced and extended zone schemes, along 
with a third, the periodic zone scheme, is given by Harrison (ref. 11). 

These different ways of picturing the wave-number distribution of the electrons 
must give identical results for the electronic properties of the system, and the choice of 
which scheme to use is dictated primarily by considerations of convenience and clarity. 
The convention adopted in this report is that integrals over wave- number space (which 
will ultimately become integrals over the Fermi surface) are considered to be over a 
single Brillouin zone. Whenever it becomes necessary to explicitly define t 1 ' region of 
integration, therefore, the reduced zone scheme shall be used. 

The most fruitful attempts to describe the transport processes in metals by means 
of the BE have started with the linearization of equation (5). The distribution function 
f(k) is expanded in a power series in the electric field, and all terms beyond the first 
power of E are dropped. The deviation from equilibrium then is written as (ref. 7) 

f(k) - f°(ej = g(k) — (8) 

9e k 

The linearization results in the following form for the BE: 

eE * v k + — v k x H • V k g(k) = J Q(k, k’)[g(k) - g(k’)]dk’ (9) 

Re 

Most solutions of this equation have started with the assumption that the scattering term 
may be simplified by means of a relaxation-time approximation. This approximation 
was used by R. G. Chambers (ref. 12), for example, who showed that there exists an 
expression for the deviation of the distribution function from equilibrium if the scatter- 
ing is governed by a relaxation time, that is, if the probability of a scattering occurring 
in a time dt is proportional to dt with the proportionality constant l/r(k). This is 
equivalent to the identification 
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( 10 ) 


jfo -_ f V] 

r(k) 

Budd (ref. 13) has shown that Chambers’ solution to the BE (which was derived originally 
to include only an electric field) needed only slight modification to include a magnetic 
field, and that the result was an exact solution to the linearized, relaxation-time BE. 

The resulting solution is in terms of integrals over the past history of the electron, and 
its exact form is not of interest here. Although the form of the solution makes it difficult 
to apply easily to the properties of real metals, applications have been made to the 
theory of transport properties of semiconductors (ref. 14), the attenuation of ultrasonic 
waves in metals (ref. 15), and the asymptotic form of the magnetoconductivity tensor for 
metals (e.g. , refs. 16 and 17). 

The introduction of a relaxation time for the scattering is primarily a device for re- 
ducing the BE (in the form of eq. (9)) to a differential equation for the distribution func- 
tion. The applied magnetic field drives a conduction electron in an orbit around the 
Fermi surface in a plane perpendicular to the direction of the applied field. Such a 
simplification of the BE has a natural set of coordinates in terms of energy, relative 
position on an orbit, and the position of the orbit on the surface. This description allows 
the contribution of an electron to the transport current to be interpreted in terms of the 
details of the Fermi surface sampled in its trajectory. Thus, there is a natural frame- 
work available to consider the effects of closed orbits, orbits that are closed but extend 
through many zones (in the extended zone scheme), or even trajectories that extend in- 
definitely through k- space, as in the case of copper. The effect of each of these types of 
trajectories may be considered as a function of the lifetime of an electron in its orbit 
before a scattering occurs. 

This approach to the transport problem was carried far by the work of Lifshitz, 
Azbel', and Kaganov (ref. 18) and Lifshitz and Peschanskii (refs. 19 and 20). These 
authors cataloged a large number of possible Fermi surface shapes in terms of their 
contributions to the transport properties in high and low magnetic field strengths. This 
catalog of qualitative descriptions still provides the working tools for the experimentalist 
who is investigating the shape of a Fermi surface from galvanomagnetic experiments. 

The preceding discussion was concerned with work which sacrificed any attempt at 
a detailed description of the scattering, the result being a problem which is mathemati- 
cally simple. Although the usefulness of the results, in terms of understanding the ef- 
fects of Fermi surface topology on transport properties, justified the relaxation-time 
approximation, this approach is not adequate to support current experimental work. The 
task of measuring Fermi surface dimensions has been taken over by experiments relying 
upon other phenomena, such as the de Haas - van Alphen effect. The galvanomagnetic ef- 
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fects, on the other hand, are more useful for investigations of the scattering of electrons 
in metals. Clearly, the relaxation-time approximation is not adequate for such investi- 
gations, and more general theories must be developed. 

One approach to a more satisfactory treatment of the scattering term in the BE was 
introduced by Sondheimer (ref. 21). He made the ad hoc assumption that the scattering 
transition probability could be written as a finite sum of products of a function of k with 
a function of k’, that is, 


N 

Q(k,k T ) = ^ b.-q.^q.^OSCej^, - e k ,) (11) 

i, 

J=1 

The matrix b- is real and symmetric, and the real functions q.(k) are arbitrary, ex- 
cept that they must be square-integrable in the Brillouin zone. In the theory of Fredholm 
integral equations, such a kernel is called degenerate, and an extensive theory of such 
integral equations has been developed (Mikhlin (ref. 22) and Courant and Hilbert 
(ref. 23)). The effect of this replacement (eq. (11)) in a Fredholm equation is that the 
integral equation is replaced by a finite set of linear equations for the matrix elements 
b^j, if the functions q^ are chosen properly. The effect of this assumption on the BE 
is to reduce equation (9) to a finite set of coupled differential equations. With this re- 
duction, Sondheimer (ref. 21) and Jones and Sondheimer (ref. 24) were able to discuss 
the resistivity and asymptotic longitudinal magnetoresistance of metals with anisotropic 
properties. Although the use of the degenerate kernel was an arbitrary construction, a 
partial justification of this form for some types of scattering has been made by Mann 
(refs. 25 to 27). A further discussion of this theory and its relation to the present work 
is made in the section DEGENERATE KERNEL THEORY. 

The approaches to the transport problem discussed here have had one feature in 
common - the linearization assumption of equation (8) or its equivalent. This linearization 
assumes that the deviation from equilibrium of the distribution function is proportional 
to the magnitude of the electric field. The starting point for the present work, on the 
other hand, is the observation by Taylor (refs. 8 and 9) that as a consequence of the 
linearization assumption equation (8) may be rewritten, displaying the electric field de- 
pendence explicitly, as 

m - I°(e k ) = -eE • A(k) ‘ ,! ° (12) 

8E k 
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Thus, the scalar function g(k) is replaced by a vector function A(k), which has the 
properties of a vector mean free path. It will be shown that the BE has an exact, closed- 
form solution for A that is valid for any strength of the applied magnetic field. This 
formalism provides a natural framework within which the transport properties of a real 
metal may be calculated. 

The fundamental assumption in the work described herein is that the independent- 
particle Boltzmann- equation description is valid for transport processes involving elec- 
trons in metals. Three basic difficulties with this assumption are the many-body aspects 
of the electron- electron interaction, the inherently nonequilibrium (in the thermody- 
namic sense) nature of the system, and the quantization of the electron system into 
Landau levels by the applied magnetic field. These difficulties may not be resolved 
rigorously from the present state of solid-state theory, but we may present points to in- 
dicate that this approach is at least a reasonable one to take. 

The questions raised by the many-body aspects of the electron- electron interactions 
in metals apply, at least in large part, to the existence of a Fermi surface in such sys- 
tems. The long-range correlations of the electrons could be considered too strong to 
allow such a sharp discontinuity to exist. It appears in reality that such is not the case. 
Work by Migdal (ref. 28) (and discussed by Abrikosov, Gorkov, and Dzyaloshinski 
(ref. 29)) shows that there is a discontinuity in the electron momentum distribution in an 
electron gas, even in the presence of interactions. Furthermore, the existence of such 
sharp experimental evidence concerning the size and shape of the Fermi surfaces in a 
large variety of metals must be considered good evidence that this construction has 
some validity. 

Whether the Boltzmann equation provides a valid description of transport processes 
in metals is a question that has received a great deal of attention. The distribution 
function for an electron in an independent-particle model is an approximation to a den- 
sity matrix and a quantum-mechanical master- equation description of a very complex, 
time-varying system. Much effort has gone into attempts to understand the role of the 
Boltzmann equation in the framework of more general formulations of nonequilibrium 
thermodynamics (e.g. , the work of Prigogine (ref. 30) and Edwards (ref. 31)). The role 
of the BE as it applies to steady- state phenomena has been extensively discussed by 
Dresden (ref. 1). Basically, it still remains an assumption that the Boltzmann equation 
is an adequate approach to steady-state processes such as electron transport in metals. 

There are a large number of experiments which observe effects in metals which are 
quantum mechanical in origin. These include the de Haas - van Alphen oscillations in 
the magnetic susceptibility, de Haas - Shubnikov oscillations in magnetoresistance 
measurements, magnetostrictive oscillations in the physical dimensions of a metallic 
crystal, and many others. These effects all arise from the quantization of electron 
orbits in the presence of the magnetic field (i. e. , the formation of the Landau levels in 
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the momentum distribution). In contrast with this list of quantum effects, there is the 
equally valid observation that, in most cases, these quantum effects are oscillations 
superimposed upon usually monotonic semiclassical effects. For example, the oscilla- 
tions which appear in the longitudinal magnetoresistance at high fields (the de Haas - 
Shubnikov effect) are structure imposed on a monotonic magneloresistance which can be 
calculated on the basis of the BE. 

Some justification of the applicability of the semiclassical method may be obtained 
from a consideration of the energies involved. Three useful energies for conduction 
electrons in a magnetic field are the Fermi energy e F , the width of the thermal broad- 
ening of the Fermi surface k-gT, and the separation of the quantized Landau levels 
Bw =KeH/mc. It is convenient to take some representative numbers for comparison; for 
Pippard's ’’standard metal" (ref. 3) these energies have the following values: 

e F = 5. 57 eV 

k B T = 3. 62x10” 4 eV(at 4. 2 K) 

Rco = 1. 16x10” 5 eV(at 0. 1 T) 

Thus, the Fermi energy is many orders of magnitude higher than either the thermal 
energy or the Landau level separation, and is the dominant energy factor. Fur- 
ther, even at 1 tesla the level separation is only comparable to the thermal broadening. 
This latter comparison is further substantiated by the observed fact, mentioned pre- 
viously, that even when quantum effects are present, there is still a substantial contri- 
bution which is of semiclassical origin. 

In summary, then, although the use of the Boltzmann equation cannot be rigorously 
justified, there is every indication that it provides a valuable tool for the description of 
conduction processes in metals. On this basis then, we shall proceed to the formalism 
for the solution of the BE for the galvanomagnetic effects in metals. 

The present work is divided into three principal segments. The following section 
takes the assumption of a vector mean-free-path for the Boltzmann equation, including 
an applied, uniform magnetic field, and traces the development of an exact, closed- 
form solution for the conductivity tensor. As indicated previously, this solution may, 
in principle, be found for all forms of the Fermi surface and the scattering model; 
elastic scattering has been assumed for simplicity of the mathematics. The next 
two sections treat the implementation of the theory, the insight that a knowledge of 
the pertinent lattice point group may bring, and a sketch of the numerical procedures 
involved in performing a calculation. Finally, results are shown for two scattering 
models. The properties of an isotropic scattering model for a free-electron gas are 
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discussed. And, then, the general implications of a degenerate kernel, as proposed by 
Sondheimer, are discussed with sample results for a spherical Fermi surface model. 


SOLUTION OF THE BOLTZMANN EQUATION 

The introduction considered the Boltzmann equation for conduction electrons in a 
metal. Steady, homogeneous electric and magnetic fields are applied to the metal, and 
the sample is assumed to be homogeneous with a uniform temperature distribution. In 
this situation, the distribution function is a function of k alone, and is written as f(k), 
or simply f^. When the scattering is elastic, the Boltzmann equation is 

^I + lv k xlJ. V k f k = f Q(k, k») (f k - f k ,)dk' (13) 

where the function Q(k, k') is the (positive) transition probability per unit time for 
scattering from the state k into the state k\ For this case of elastic scattering, 

Q(k, k’) is symmetric (i. e. , Q(k, k') = Q(k',k)) and is proportional to the energy- 
conserving 6(e k - e k ,). 

The objective of the present work, as discussed in the introduction, is a solution of 
equation (13) that does not depend upon the precise nature of the Fermi surface or the 
details of the scattering mechanisms. In this section a closed-form expression for the 
conductivity tensor is found. (An earlier version of the formalism described in this sec- 
tion has appeared as ref. 32. ) 

The equilibrium distribution function for the electrons is the Fermi function f°(e k ). 
An applied homogeneous magnetic field does no work on the electronic system, and the 
effect of an electric field is to perturb the distribution function to a new function f k . 

The deviation from equilibrium is assumed to be an expansion in powers of E (refs. 8 
and 9), that is, 


f k " = f(1) ( E ) + f (2) (E 2 ) + . . . (14) 

The linearization of the BE is obtained from the assumption that only the first term 
f^(E) of equation (14) is important. This assumption results in a transport current 
which is also proportional to E, the condition of Ohm's law. The explicit form of the 
linearization is (refs. 8 and 9) 
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The linearized Boltzmann equation for the function A(k) is 


v k ‘ — < v k X H ' = f k 't A M ~ A(k’)]dk f (16) 

he • / 

The basic quantity describing the scattering is the positive, definite, symmetric 
kernel Q(k,k’). It will be convenient for the later discussion to define an integral oper- 
ator Q such that 


Qx(k) s J' Q(k, k^x^^dk* (17) 

and a scalar function 

q(k) h J Q(k,k’)dk’ (18) 

With the introduction of these quantities the Boltzmann equation may be written as 

|(q - Q) + ~ (v R x H • V k )j A(k) = v k (19) 

The essence of the method of Taylor (refs. 8 and 9) is the use of the operator q - Q 
to treat the scattering. No assumptions about the topology of the Fermi surface or the 
details of the scattering processes are needed for the solution of the BE. It has been 
shown by Taylor (refs. 8 and 9) that the operator q - Q possesses an inverse, as long 
as the scattering kernel is not a 6-function of k and k', and if the constant function is 
excluded from the range of the operator. Both conditions can be shown to be satisfied 
for all problems of interest. We define the differential- integral operator 

T = — (q - O" V k * h x V k ) (20) 

Tl 

A «*■ 

where h is the unit vector in the direction of the magnetic field (i. e. , H = Hh). Then, 
operating on the left side of equation (19) by (q - Q)~* and defining the free-electron 
cyclotron frequency 



result in 


( 21 ) 


me 


(1 - io>T)A = (q - Q) _1 v k 

Equation (21) represents the Boltzmann equation for the function A(k). If a solution 
can be found for A, then the transport properties of the system may be found from equa- 
tion (15) for the distribution function. For a system with no applied magnetic field (on - 0 
in eq. (21)) it was shown by Taylor (refs. 8 and 9) that there exists a solution for A in 
terms of eigenfunctions of the operator q~*Q. For a nonzero magnetic field, equa- 
tion (21) can be solved by forming the operator which is the inverse to 1 - iwT and op- 
erating on the left side. The inverse operator (1 - iwT) - * exists since T, although not 
Hermitian, possesses a real eigenvalue spectrum (see appendix A). Hence, 1 - iu>T 
cannot be singular (ref. 33) and the inverse exists. The formal solution for A then is 

A(k) = (1 - iu)T) _1 (q - Q) -1 v k (22) 

The transport properties of the system, and in particular the conductivity tensor, 
may be determined by substitution of equation (22) into the distribution function (eq. (15)). 
The conductivity tensor, the quantity that describes the galvanomagnetic properties of 
the metal, is 


a = 



(1 - ia>T) 


-1 (q - Q) 



dr 


dk 


(23) 


Equation (22) represents, symbolically, an exact, closed-form solution to the lin- 
earized Boltzmann equation. Thus, the conductivity tensor (eq. (23)) is an exact result, 
to within the linearization approximation, and is valid for any strength of the magnetic 
field. This generality should be contrasted with previous theories which were formu- 
lated to describe either low-field or high-field behavior. No assumption about the nature 
of the Fermi surface has been required, and the restriction to elastic scattering may be 
removed by the use of an inelastic scattering operator P (refs. 8 and 9) rather than the 
Q of equation (17). 

It is apparent that equations (22) and (23) are not in a form suitable for computation. 
To get a useful form for the conductivity, the quantity of interest here, we must con- 
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sider the meaning of the two inverse operators in equation (22). It is shown in appendix A 
that the operator T, although not Hermitian, does possess a spectrum of real eigen- 
values and a corresponding complete set of (in general, complex) eigenfunctions. These 
eigenfunctions are referred to as the "magnetic scattering eigenfunctions" and form a 
complete set for functions in k-space. These eigenvalues and eigenfunctions are defined 
by the equation 


Tb z (k) = T z b z (k) 


or, explicitly, 


— ( y k • h x V k )b z (k) = r z (q - Q)b z (k) (24) 

li 

It follows from equation (24) that Tb* = -i^b* so that the b* are the eigenfunctions of 
T corresponding to negative 

In both equations (22) and (23), the operator T acts upon the zero-field mean free 
path (refs. 8 and 9) (q - Q)” 1 v^. Since the {b z } form a complete set, the following ex- 
pansion can be written: 


(q - <D'\ - 


( 25 ) 


l 


th 


where v is the p Cartesian component of the velocity vector v, . The expansion 

IX K 

coefficients /3^ z , although they do not depend upon k, do depend upon the vector h, the 
direction of the magnetic field, from the definition of T. Note that, since the left side 
of equation (25) is a real function, for every term there must also appear the 

corresponding /3* z b*, guaranteeing the reality of the sum. 

With this expansion for the zero-field mean free path, the conductivity can be 
written as 


a 


v\i 




(26) 


The presence of the energy derivative of the Fermi function in the integral of equa- 
tion (26) leads to some discussion of the k-dependence of functions in this theory. Be- 
cause of the limit 
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= - 6 ( e k - e F> 


Lim 

T-0^ 


only states at the Fermi surface can contribute to the integral in equation (26). The im- 
portance of this energy 6-function may be seen by consideration of a variable transfor- 
mation (ref. 7) from the volume element dk to an equivalent volume element 


(27 ry 


3 v de k dS 

d^k = — *_ 

(2tt) 3 R l v kl 


(27) 


where dS is the element of area on the energy shell corresponding to This sep- 
aration of variables into energy and surface variables, combined with the 6(6^ - ep) of 
equation (26), implies that only functions evaluated on the Fermi surface need be con- 
sidered. To simplify the notation in integrals, it is convenient to redefine the volume as 

dk s de^ df2^ 


where dfi^ stands for the surface integral and the multiplying factors in equation (27). 
After the integrations have been performed over the energy 6-function, all that remains 
is the integral over dfi^, or simply dfi, which represents an integral over the Fermi 
surface. 

It is important to make the connection between the preceding discussion and the re- 
striction to elastic scattering. The transition probability Q(k, k') is defined for every 
k and k' in wave-vector space. The assumption of elastic scattering restricts atten- 
tion, for a given k, to only those k’ that lie on the same energy shell. Since every k 
is considered, the factoring of k- space into a nested set of constant- energy surfaces is 
a natural view. There is no further restriction implied by the condition of elastic scat- 
tering. The functions of interest, such as Q(k, k'), q(k), and A(k), are defined off the 
Fermi energy surface. The discussion of the preceding paragraph is a low-temperature 
approximation; that is, if these functions are slowly varying functions of energy near 
the Fermi energy, then all k-vectors may be restricted to the Fermi surface. 

After equation (26) the problem becomes the calculation of the field- orientation- 
dependent eigenfunctions b^(k) and the expansion coefficients )3 It is apparent from 
equation (24) that the solution for the { b^ } involves an integrodifferential equation. For- 
tunately, the prospect is not as dismal as it first appears, for in equation (24) the inte- 
gral operator Q acts upon the function b^. Thus, if the {b^} themselves could be ex- 
panded in terms of the eigenfunctions of q _ *Q, the problem could be simplified. This 
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simplification turns out to be quite substantial and illuminating. 

Taylor, in his treatment of the magnetic-field-free Boltzmann equation (refs. 8 
and 9), showed that eigenfunctions of the non-Hermitian operator q” Q, 

(q _1 Q)a n (k) = <* n a n (k) (28) 

could be defined. Appendix A includes a discussion of the properties of these scattering 
eigenfunctions, with their existence and completeness properties. Among their useful 
properties is the fact that they may be chosen to be orthonormal, so that 

f a r 0«)q(k)a n (k)df2 = 6 rn (29) 

These scattering eigenfunctions are real, and depend only upon the nature of the Fermi 
surface and the electronic scattering processes, independent of the magnetic field vec- 
tor. 

The magnetic scattering eigenfunctions may be expanded in the a n (k) as follows: 

(30) 

n 

where the complex expansion coefficients y^ n depend implicitly upon the field direction 
h. With the aid of this expansion and the orthogonality condition of equation (29), the 
eigenvalue equation for the becomes 

-imh • £ y ln /a r (k)[v k x V k a n (k)]df 2 = Kr^l - a r )Y lr (31) 

n 

Although this equation appears formidable, its true nature can be seen by a definition of 
a new set of coefficients 



and the vector quantity 

C r „ * - Y [a - %>(1 - “r>]" 1/2 / X V k a n«> < 33 > 
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Thus, equation (31) has the form 


fc-Z C rn r !„ = V;r < 34 > 

n 

This equation represents a set of simultaneous equations, in effect a secular equa- 
tion, for the The important feature of this equation is the fact that the only depend- 
ence upon the magnetic field is through the presence of h. The matrix C is determined 
completely from the field- independent scattering eigenfunctions a n (k). 

It can be seen by an integration by parts of equation (33) that the C matrix is 
Hermitian and, in fact, purely imaginary. This implies that the eigenvalues must 

be real. Also, since the trace of C must vanish for any field direction, the must 

occur as positive-negative pairs; that is, the sum of the must also vanish. 

The combination of equations (26) and (30) gives as an intermediate form for the 
conductivity 


where 







^JjjTlX^vn 
1 - iCOTj 




/ 




(35) 


(36) 


The final reduction of the conductivity is motivated by the realization that equa- 
tion (25), the expansion of (q - Q) _1 v in terms of the b 7 , is equivalent to an expansion 
of the same quantity in terms of the scattering eigenfunctions a n (k). This, in fact, was 
the procedure used originally (refs. 8 and 9) in the solution of the BE for the zero- 
magnetic-field conductivity. This double expansion results in 


<« - =1 V*® =11 < 37 > 

l In 

The a n are eigenfunctions of q”^Q, so operating on both sides of equation (37) with 
q - Q and using the orthogonality of the a n result in 
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g M r = Y, P m Y lr (1 " 

l 

with the g defined as in equation (36). Since C _ is Hermitian, the r ; form an 
j£T rii h /o t-t 

orthonormal set of vectors. The replacement of y^ r ( 1 - Q! r ) ' by r^ r , as in equa- 
tion (32), and the use of the orthogonality condition result in an expression for ®\lI 



E* 

Vlr 


This result allows the final simplification of the expression for the conductivity. From 
equation (35), then, 


cr 


vii 


g 2 \ ' 1 \ % g pr r £r r £n g im 


(38) 


or 


vii 


6 1 _ iwr 2 g Mr 7 fr^n g i/n 

-S— ' 1 7 r, n 


(39) 


It is not immediately apparent that either form for the conductivity is closely akin 
to the usual results (ref. 3) presented on the basis of phenomenological arguments. Any 
detailed comparison would involve some specification of a model for the Fermi surface 
and would not have the general validity of the preceding discussion. The desired com- 
parison may be helped by the recognition that the magnetic scattering eigenvalues 
occur in positive- negative pairs. In fact, as may be verified from equations (33) 
and (34), if r^ n is an eigenvector corresponding to an eigenvalue r^, then r* n is also 
an eigenvector with the associated eigenvalue -t^. Thus, if the sum over l in equa- 
tion (39) is restricted to only include the positive values of t^, the final result is 
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SYMMETRY PROPERTIES OF THE SCATTERING EIGENFUNCTIONS 

It was shown in the preceding section that a formal, closed-form solution to the 
linearized Boltzmann equation exists for all magnetic field strengths and orientations. 

It was also demonstrated that this solution could be obtained by using two eigenfunction 
expansions. For any application of these results to a real metal these eigenvalue prob- 
lems must be formulated in terms of the electronic scattering processes and band 
structure. This places a limit upon the number of general observations which may be 
drawn from the formalism. Fortunately, a large amount of information about the scat- 
tering eigenfunctions (referring now to both eigenvalue problems) and the methods of 
computation may be obtained from the symmetry properties of the system. In this sec- 
tion the application of group theory to determining the properties of the scattering eigen- 
functions is discussed. 

The difference between a free-electron gas and the system of conduction electrons 
in an independent- particle model for a real metal is the anisotropy forced upon the con- 
duction electrons by the lattice of the ion cores in the metal. For many metals the most 
striking example of this anisotropy is the shape of the Fermi surface, the constant- 
energy surface which, at absolute zero, defines the boundary between filled and unfilled 
electronic states in k-space. For some metals, such as aluminum, the Fermi surface 
is very complex and is composed of several separate sheets; for copper, the Fermi 
surface is a single sheet but is distorted from the free-electron sphere sufficiently to 
touch the boundary of the Brillouin zone. This contact provides a path through k-space 
for electrons, giving rise to the concept of an "open" orbit for conduction. The common 
ground between these different cases is the influence the symmetry of the lattice has 
upon the electrons. The lattice has the symmetry properties of a specified point group, 
and the Fermi surface also has the properties of the same point group, with respect to 
rotations and rotatory inversions in reciprocal space. The properties of the conduction 
electrons must also be governed by this point group, and the scattering eigenfunctions 
may be considered from this viewpoint. 

Consider a Fermi surface whose symmetry properties are those of a point group S?. 
This group consists of some number g of rotations and possibly rotatory inversions, 
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each operation symbolized by the operator St- This group has a number of irreducible 
representations which give a set of finite matrices to represent the operations. If the 
irreducible representation of & has the dimension of n^, then the matrices for 
this irreducible representation are square and have dimension n^, and the matrix ele- 
ments are represented by the symbol D^(0?). The effect of a group operation acting 
upon a function of k is defined by 

#f(k) = f(R -1 k) (40) 

where St is symbolic for the abstract group operator and R is the 3x3 matrix repre- 
sentation of St which acts upon the triplet of numbers (k^ k^, k z ). 

Since the scattering in the linearized approximation is independent of the magnetic 
field, the Q(k, k T ) must reflect the symmetry of the lattice point group. The only as- 
sumption about the scattering kernel which must be made is a weak restriction, particu- 
larly in view of the restrictions other authors (ref. 24) have made. Consider the scat- 
tering from a state k to (or from) a state k'. If there are no preferred directions in 
space imposed by some external source, then this scattering process should be identical 
to the scattering from Rk to Rk', where k and Rk are equivalent points. The sense 
of equivalence here is simply that the point Rk is obtained from the point k through 
one of the group operations St, Thus, for the group of order g, a general point on the 
Fermi surface will have g distinct images under the group operations, including the 
identity operation. This equivalence assumption about the scattering kernel may be for- 
mally written as 


Q(k, k*) = Q(Rk, Rk’) (41) 

for all k and k' and all St^W. This assumption does not require that Q(k, k T ) be de- 
generate or that it only depend upon k - k' and, in fact, is a very weak restriction upon 
the form of the scattering kernel. The fact is, the use of this assumption and the power- 
ful tools of group representation theory allow interpretation and understanding of the 
scattering eigenvalue problem (eq. (28)). 

For the purposes of the problems considered here, the most important result from 
group representation theory is the concept of a basis function of an irreducible repre- 

1L 

sentation. The formal definition of a basis function is as follows: If the cr n irreduc- 
ible representation of a group $ has the dimension n . then there exist n, basis 
( a \ a a 

functions ’ such that, for any & in 57, 

^\ a) (k) = Yj ^° /) (k)D^ ) (R) (42) 
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To make the connection between the scattering eigenvalue problem and the basis 
function the following theorem is important: 

Theorem I : The integral operator Q of equation (22) commutes with all St of & , that 
is, 


#(Qf) = Q(5K) 


(43) 


Proof : By the definition of Q, 

#(Qf) = St f Q(k,k’)f(k’)dTC’ 

= J' Q(R - ^ k, k')f(k')d£2' from equation (40) 

= y Q(k, Rk^ffk^dn' from equation (41) 

but f(k') = f(R~*Rk') = £Pf(Rk') and the energy of Rk' is the same as the energy of k’ 
due to the invariance of the Fermi surface. Thus, 

<*(Qf) = /q£, Rk')<*f(Rk')df2’ 

= Q(s*f) 

All k vectors have been taken to be vectors to the Fermi surface, to be consistent with 
the discussion of the preceding section. All integrals are taken over the Fermi surface. 
This assumption is not as restrictive as it appears, since the energy integral implies 
an integral over all the energy shells. Each constant- energy shell has the same sym- 
metry properties as the Fermi surface, and, thus, the generalization is immediate. 

Theorem I implies that there is a connection between the scattering eigenfunctions 
of equation (28) and the set of basis functions of the group representations. The estab- 
lishment of this connection, along with further implications, is the object of this section. 

The scattering eigenvalue problem (eq. (28) or, in practice, eq. (Al)) is in reality a 
Fredholm integral equation, a homogeneous equation of the second kind (ref. 22). It is 
a result of a theorem in the theory of integral equations known as Mercer's theorem 
that, since Q(k, k’) is a positive definite, continuous, symmetric kernel, there exists 
a biharmonic expansion for Q in terms of its eigenfunctions (ref. 23). The translation of 
this theorem requires a little groundwork. 

The spherical harmonic functions Y* 11 (6 , <p) form a complete orthonormal set of 
functions for functions defined on the unit sphere. Another property of the spherical 
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harmonics, perhaps not so widely recognized, is that they are basis functions of the ro- 
tation group. The Y™ obey an equation analogous to equation (42) (ref. 34) with Z 
labeling the irreducible representations and m denoting the partners of the basis. 

Since the Z u irreducible representation of the rotation group has dimension 21 + 1, 
there are 21 + 1 partners in a basis set, corresponding to the 21 + 1 values of m for 
a given Z in a spherical harmonic. 

Now consider a kernel Q(k, k'), where k and k’ are vectors to the unit sphere. 

If this kernel is restricted to being spherically symmetric (i. e. , the function Q may 
depend only upon the difference |k - k’ f), then there exists an expansion of Q in terms 
of the spherical harmonics 

Q(k,k’) = Yj C l ^ m (k)Y z m *(k') (44) 

l~ 0 m 

Note that the existence of this expansion does not imply that the kernel must be degen- 
erate, in which case the sum would have to be finite. The implication in equation (44) is 
that the function Q is well approximated by the sum on the right side if a sufficiently 
large number of Z-values are included. 

The purpose is to find an expansion, similar to equation (44), for a scattering kernel 
which is restricted to having the symmetry of a finite point group, as specified by equa- 
tion (41). Obviously, the functions analogous to the spherical harmonics (i. e. , a com- 
plete, orthonormal set of basis functions of the point group S?) are required. 

The existence and methods of obtaining such functions, known generally as 
symmetry-adapted functions, are sufficiently detailed in the literature (ref. 35) not to 
require an exhaustive treatment here. Given a group ^ (meaning that we have all the 
matrices of all the irreducible representations of S7), consider one of the irreducible 
representations, labeled by a, which is of dimension n . Since g is assumed to be 
a point group, it is a subgroup of the full rotation group, including inversions, if neces- 
sary. The representations of the rotation group are, in general, reducible with respect 
to <<?. This means that the Z^ representation of the rotation group may contain more 
than one representation of <& . If the irreducible representation is contained in the Z** 1 
representation, then there are n ff linear combinations of Y™, summed over m, that 
are basis functions for the a representation of The set of all such linear com- 
binations, for all Z and all a of S7, form a complete orthonormal set of functions. 

The symmetry- adapted functions for the full cubic group have long been known 
and tabulated as the Kubic harmonic functions (ref. 36). In principle, the same con- 
structions may be performed for each crystal group (in particular, see Altmann and 
Cracknell (ref. 37)) and all such tabulations could be kept as ready reference. How- 
ever, all that is required is the existence of such a complete orthonormal set of func- 
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tions X^ r indexed by the angular momentum of the contributing spherical harmonics l, 
the label of the irreducible representation of & j, and the column (or the component) of 
the representation r. 

Consider a kernel function Q(k, k'), where k and k' are vectors to the Fermi 
surface, which is to be invariant under the operations of a point group, as discussed 
earlier in this section. Then there exists an expansion for Q, analogous to equation (44), 
in terms of the symmetry-adapted functions. This expansion is the content of theorem n, 
as follows: 

Theorem II : Given a function Q(k, k') which is invariant under all operations dt of a 
group S? (i. e. , Q(k, k’) = Q(Rk, Rk') for all in 9), then there exists the expansion 

Q(k,k’)= Yj (45) 

i,v, i r 

The proof of this theorem is straightforward, but not instructive for the purposes of this 
report, and so it is given in appendix B. The important fact to note is that equation (45) 
does not represent a treatment of the scattering by a degenerate kernel, since the sums 
of l and V are infinite. Equation (45) does however, indicate, simplifications which 
arise if the kernel is degenerate. Further evidence of this point appears in the section 
DEGENERATE KERNEL THEORY. 

The groundwork is now laid to make the connection with theorem I. It is apparent 
that the eigenfunctions of the scattering operator must be related to symmetry- adapted 
functions. This connection is provided by the content of theorem HI: 

The orem HI : A function yf(k) is a basis function for the group <&; that is, it satisfies 
the equation 


Sty] (k) = Y yj' , (k)D ] rrT (R) (46) 

r' 

if and only if 

y 3 r (k) (47) 

l 

The proof of this theorem is also given in appendix B, since it is merely an exercise in 
group representation theory. Equation (46) is the same as equation (42), although the 
notation has been changed to make the following discussion smoother. 


L 
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Now that the connection between the symmetry-adapted functions and the scattering 
kernel has been established, the scattering eigenvalue problem may be reformulated. 
Since any kernel Q may be expanded as in equation (45), the modified scattering kernel 
of equation (Al) of appendix A clearly has such an expansion. It is shown in appendix B 
that the function q(k) is a scalar with respect to the group operations, so that the modi- 
fied scattering kernel has the same symmetry properties as the original kernel Q(k, k T ). 
Formally then, 


-f=== £ 4.£*r«^y> (48) 

f q(k)q(k') l,V, j r 

The relation between the basis functions yF and the symmetry-adapted functions X^ r 
(eq. (47)) is used to derive a secular equation for the scattering eigenvalues a, as in 
equation (Al), 


X/ ~ a 6 u) b r ~ 0 ( 49 ) 

v 


A consequence of theorem n, as shown in appendix B for real, symmetric kernels, is 
that the matrix of expansion coefficients A must be Hermitian. Thus, the scattering 
eigenvalues, the roots of equation (49), must be real, confirming the results of appen- 
dix A; and the vectors of the expansion coefficients, the b^, must either be automati- 
cally orthogonal (for distinct eigenvalues), or else be constructed to be orthogonal (for 
degenerate eigenvectors). Since the scattering eigenfunctions, labeled by a*" according 
to the notation of equations (47) and (28), are related to the eigenfunctions of equation (48) 
through the scalar function Vq(k), the transformation properties of the two functions 
will be the same. 

The meaning of the notation and the scattering eigenvalues can be seen by looking at 
the structure of the matrix A of equation (48). This coefficient matrix connects only 
symmetry functions which belong to the same irreducible representation of the group. 

In this sense, A is a diagonal block matrix of block dimension g and each block, in- 
dividually Hermitian, is indexed by its irreducible representation. Each block will be 
of infinite dimension if the kernel is not degenerate and its eigenvalues are the scattering 
eigenvalues. Each scattering eigenvalue must carry two indices, one corresponding to 
the irreducible representation and one corresponding to an ordinal number within its 
block. There is one complication which is partially hidden in this indexing. If the a 
irreducible representation of S? has the dimension then each eigenvalue corre- 


24 



sponding to the block a must be n^-fold degenerate, corresponding to the n ff partners 
of the set of basis functions. 

The possibility of providing an indexing and a characterization of the symmetry 
properties of the scattering eigenfunctions allows an insight into the meaning and struc- 
ture of the formal equations derived in the preceding section. The fact that it also pro- 
vides an aid in the numerical treatment of the problem will be considered in the next 
section. For the time being, consider two major steps in the formalism, the formation 
of the C-matrix, (eq. (33)) and the definition of the g yn factors of equation (36). Use of 
the tools of group representation theory allows some strong statements to be made con- 

f 

cerning these equations. 

First, consider the definition of the g yn , but with the indexing of the scattering 
* eigenfunctions indicated explicitly 

S vl = ~ f v^a*"’ n (k)dS2 (50) 

The indices r and n on a refer to the column of the representation and the ordinal 
number of the eigenvalue, respectively. The index l on the g is meant to include 
schematically all possible indices from the a. For convenience, let us consider the 
full cubic group of order 48 with 10 distinct irreducible representations, having 
dimensions from 1 to 3. The velocity vector v^. will have specified transformation 
properties for this group; in fact, v^ transforms according to the r l5 irreducible 
representation (ref. 38). The selection rule for the integral in equation (50) may be ob- 
tained from one of two approaches. Either the integral is formally evaluated by sum- 
ming over the contributions of the 48 equivalent zones of the Fermi surface, or else the 
content of the Kronecker product ^(v)®^(aj) is examined for a scalar part. The 
script is used to symbolize the irreducible representation of the indicated quantities. 
The former approach is used in the following section, but the latter is more useful for 
the investigation of the selection rules. 

f The appropriate guide to use is the result (ref. 39) that an integral of the type shown 

in equation (50) will vanish unless there is some part of the integrand which is a scalar 
under the group operations. Since, for O^, v^ is known to transform as Tjg, all that 
is necessary is to scan a table of Kronecker products of Tjg with all irreducible rep- 
resentations of to discover the fact that only r l5 contains the identity Tj. 

Even further, the integral will vanish for a scattering eigenfunction which is unless 

the column index of the function is the same as the column index of the velocity (i. e. , the 
component of v^). Thus, 



af> n (k)dJ2 6. 


],9{vfrv 


(51) 
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By the same token, the selection rules for the formation of the C-matrix (eq. (36)) 
can by obtained by looking at a triple Kronecker product of representations. This should 
not be striking since C is really the matrix of the matrix elements of a crystal angular 
momentum operator v k x V taken with respect to the scattering eigenfunctions. This is 
a familiar problem in nuclear physics, and its solution is the set of Racah coefficients 
for the scalar content of the triple product. 

In the full cubic group the angular momentum transforms as the Tjg irreducible 
representation and thus has even parity under inversions. A table of the triple product 
of Tjg with all the irreducible representations will show which products contain the 
identity representation. Conversely, for the products which do not contain the identity, 
the matrix element of C must vanish. Even without such a product table, the angular 
momentum cannot have matrix elements between states of different parity, since the 
product must then have odd parity under inversion. Thus, the only contribution to the 
gross transport properties arises from the part of the scattering kernel which has odd 
parity, because of the nature of the g^. The qualification implied here is necessary, 
since the effective anisotropic relaxation time q *(k) is determined from the scalar part 
of Q(k, k') and determines the scale for the system. 

The importance of the symmetry properties in determining the transport properties 
of electrons in a metal should be apparent at this point. The representation theory of 
finite groups provides powerful tools for the understanding of the formalism discussed 
in the preceding section. The validity of an even stronger statement, that the use of 
group theory is essential for the numerical description of the galvanomagnetic properties 
of a real metal, is discussed in the next section. 


STRUCTURE OF THE CALCULATIONS 

The desired product of this investigation is a computational procedure capable of 
giving the galvanomagnetic properties of any metal. Equations (38) or (39) describe the 
conductivity of a metal, on the basis of an exact solution of the Boltzmann equation. The 
implementation of this formalism to a real metal involves a number of computational 
problems, which are considered in this section. 

There are two primary areas for comparison with experimental results, the Fermi 
surface construction and the scattering model. When a galvanomagnetic experiment is 
used as a topological probe of a Fermi surface, such as that of copper (ref. 40), the 
details of the scattering model are not as important as an accurate construction of the 
Fermi surface and velocity vectors. On the other hand, as mentioned in the introduc- 
tion, galvanomagnetic experiments have been supplanted, for the most part, by other 
types of experiments for topological investigations. The ability to investigate the effect 
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of different scattering models upon the galvanomagnetic properties of metals becomes a 
more important role. The input information to a calculation consists of a construction 
of the Fermi surface £p(k), the Fermi velocity Vp(k), and the scattering model 
Q(k,k’). 

The calculation in which equation (38) is evaluated may be divided into two distinct 
parts. It should be apparent from the discussion of the section SOLUTION OF THE 
BOLTZMAN EQUATION that the magnetic field plays no essential role in the theory until 
the determination of the magnetic scattering eigenvalues and eigenfunctions and b^ 

} is desired. The computations then may be classified as field-free or field- dependent. 

In the remainder of this section, the numerical procedures used to implement the theory 
are described in sequence. 

The computation that requires the most time and accuracy is the solution of equa- 
tion (38) for the scattering eigenvalues and eigenfunctions. Even when the scattering is 
limited to Fermi surface states, as assumed herein, the numerical problem is imposing. 
Equation (38), or, more properly, equation (Al), is a homogeneous Fredholm integral 
equation. Standard iterative procedures exist for the solution of such equations 
(ref. 22). The recognition that standard techniques are available for a solution does not 
reduce the magnitude of the numerical problem. To cover the Fermi surface with a 
mesh having a 5° separation between points would require approximately 2500 points. A 

complete knowledge of the scattering kernel Q(k, k r ), with no simplifications from the 

2 

symmetry, would then require (2500) values to be stored. It is obvious that this is a 
number larger than the fast-access storage of any computer now in use. Fortunately, 
the properties of the covering point group allow drastic simplifications in the problem. 

Let us consider, where possible, the specific example of the full cubic group of 
ordre 48. Since there are 48 elements in the group, only one forty- eighth of the Fermi 
surface is necessary to completely determine the entire surface. By the same token, 
a set of basis functions for an irreducible representation of need be specified only 
in one forty- eighth of the surface in order to compute the functions over the entire sur- 
face. It was shown previously that the eigenfunctions of equation (Al) transform as 
f basis functions of the irreducible representations of the reciprocal lattice point group. 

The degeneracy of an eigenfunction is the same as the dimension of the pertinent irre- 
1 ducible representation. This prior knowledge of the symmetry properties of the scat- 

tering eigenfunctions allows the scattering eigenvalue problem to be reduced to a prob- 
lem over one forty- eighth of the Fermi surface. Consider the eigenfunction of equa- 
i — th 

tion (Al) Yp(k), which transforms as the 1 n column of the r irreducible representation 
°f O^. The matrices of r have the dimension n^ so that the eigenvalue correspond- 
ing to this eigenfunction is n^-fold degenerate. Rather than iterating equation (1) over 
the entire Fermi surface to obtain each eigenfunction, it is only necessary to note that 
Yp is the solution to the reduced scattering problem 
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(52) 


<*Y X r (k) jf 

j=l 

The integral in equation (52) is over all k’ within one forty-eighth of the Fermi surface, 
the irreducible zone (IZ); and the kernel ^J\(k,k r ) is an n r x n^ matrix defined as 

R 

Equation (52) is an integral equation that is equivalent to equation (Al), with k and k' 
restricted to lie within the IZ. Interation of equation (52) gives the scattering eigen- 
values and eigenfunctions. 

The savings in time and storage locations that are allowed by this reformulation of 
the problem should be apparent. The convenient choice of 66 points in the IZ, equivalent 
to over 3000 points on the entire surface, entails the storage of x (66)^ values of 
■S.j rather than (3000) ^ values of Q. Since the time to perform one iteration is pro- 
portional to the square of the number of mesh points, this reduction to the IZ also pro- 
vides a substantial savings in computer time. 

The method of solution chosen is known as Kellogg’s method (ref. 22). An initial 
choice for the set of functions ■ Yp is iterated by means of the integral in equation (52). 

The eigenvalue is approximated by the ratio of the norms of successive iterates. Once 
satisfactory convergence to an eigenvalue has been obtained, a process which takes rel- 
atively few iterations, the content of the kernel due to this eigenvalue (ref. 22) is sub- 
tracted from the kernel. The iterative process then is restarted. One feature of this 
method is the use of a pseudorandom number generator to provide the initial set of func- 
tions. This choice eliminates the possibility of a trial function that is accidentally or- ^ 

thogonal to the eigenfunctions of a reduced kernel. •» 

Once the scattering eigenvalues and eigenfunctions have been obtained, the next step 
is the calculation of the C matrix (more properly, the set of three matrices which have 
been symbolized by the vector notation for C) defined by equation (33). As discussed in 
the preceding section, the selection rules for the formation of the matrix elements of C 
are obtainable from the content of the triple Kronecker product of the irreducible rep- 
resentations participating. One aspect of this observation was that equation (52) need 
only be solved for the eigenfunctions which have odd parity under inversion, since only 
these functions can give any contribution to the transport properties. 


Q(k, Rk’) 
|/q(k)q(k’) 


28 



The only numerical problem involved in the construction of the C matrices (still a 
field-free problem) lies in the formation of the gradients of the scattering eigenfunctions. 
Numerical differentiation by some type of finite-difference method is a very inaccurate 
process on a computer. To avoid this inaccuracy, the scattering eigenfunctions are ap- 
proximated by the formation of interpolation polynomials formed from the appropriate 
Kubic harmonics (the symmetry- adapted functions for the cubic group). Since the inter- 
polation polynomials may be differentiated analytically, the finite-difference method is 
not necessary. Thus, once the factors of the integrand of equation (33) are known, it 
only remains to either perform the integrals numerically for all possible products al- 
lowed by the selection rules, or else use a table of Racah coefficients for the point group 
to form the scalar parts of the integrands, with the number of numerical integrations 
being correspondingly reduced. 

The preceding discussion covers the field-free portion of the calculations. The 
numerical procedures involved are the most time- and storage- consuming of the entire 
problem. The important point to remember is that these calculations need be done only 
once for a given Fermi surface and scattering model. Once a Fermi surface construc- 
tion is made and a scattering model formulated, the scattering eigenvalue problem can 
be solved and the C matrix and the g^ n factors calculated, and a very large amount 
of information is available for use. That such a division of labor exists is clearly a re- 
sult of the linearization procedure used to simplify the Boltzmann equation. In this ap- 
proximation, the scattering term is independent of the electric and magnetic fields. To 
proceed to the next order in the Boltzmann equation would certainly introduce substantial 
difficulties. 

The remainder of the calculation, the field-dependent part, presents no essential 
problem and requires no major computational effort. Once a field direction h had been 

-A. 

chosen, the scalar product of h with the C matrix results in a purely imaginary 
Hermitian matrix which must be diagonalized. For convenience this resultant matrix is 
treated as a real, antisymmetric matrix; the imaginary parts of the resulting eigen- 
values are the (real) eigenvalues of h • C. The use of a real matrix is prompted by the 
existence of excellent methods to obtain the eigenvalues and eigenvectors. These pro- 
cedures have not as yet been extended to the case of complex matrices. The eigenvalues 

A 

are obtained by a QR transformation of h • C (refs. 41 to 43). The eigenvectors of the 
matrix are calculated from the Wielandt inverse iteration, as described by Wilkinson 
(ref. 43). These methods provide fast and accurate computations of the eigenvalues and 
eigenvectors of a general, real matrix. 

The remaining portion of the field-dependent calculation is the final evaluation of 
the conductivity tensor, and its inversion to obtain the resistivity, as a function of mag- 
netic field strength. 
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A.'i E/A'-' 'Li. 5 -Him C-i- =E3.V SURFACE A\D ISOTROPIC SCATTERING 

A csefcl crci.-_-.ci e i~ which ail caicuiarions can be performed explicitly is the prob- 
lem of isotropic scattering in a free-electron model. For isotropic (spherically sym- 
metric; scattering, the kernel Q(k, k') has the expansion 


Q(k,k’) = S y r(k) Y ”*(k’)6(e k - e k .) (54) 

Z m 

where the Y^ 1 are the normalized spherical harmonic functions. 

The integrals over k- space may be transformed into integrals over the angles 0 
and cp by the following considerations. The original definition of the volume element 
was 


(2tr) 3 

2 2 

From a transformation to spherical coordinates and the fact that e = R k /2m, 

s = _V_ mk d£ ds 
(2ir) 3 !i 2 

where dS = sin 6 dQ d<p. Since each k is restricted to the Fermi surface, the energy 
integration may be performed. Thus, 

dk — — dS = dS2 (55) 

(2/r) 3 R 2 

The scattering eigenvalue problem is solved with the aid of the orthogonality condi- 
tion of the spherical harmonics, 

Ar’< S )Y“'(S)dS = (56) 

The integration over the energy that is implicit in the operation QY™ may be done 
immediately, from the preceding discussion, and thus, 
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V 




Q(k, k* 


)Yj n (k')dfi' = 


mk 


(2tt) 3 R 2 


- 


where the point k has coordinates (k_,0, </>). Then, also 


q(k) = -L — C 

(2tt) 3 R 2 


0 


Thus, the scattering eigenvalue problem has the immediate solution 


QY z m = ^q(k)Y™(k) (57) 

with ctj = C 7 /C n . Each eigenvalue a , is (21 + l)-fold degenerate, corresponding to the 
21 + 1 functions |Y^ j for each l. The i irreducible representation of the rotation 
group has matrices of dimension 21 + 1 and the Y™ are the basis functions for the 
irreducible representation. The degeneracy of the scattering eigenvalues is due to the 
rotational symmetry of the scattering kernel, in the same manner as indicated for the 
point groups in the section SYMMETRY PROPERTIES OF THE SCATTERING EIGEN- 
FUNCTIONS. 

The usefulness of this example lies primarily in the fact that the magnetic scattering 
eigenvalue problem can be solved analytically. The definition of the magnetic scattering 
problem was, from equation (24) 

— (v k • h x V)b(k) = T (q - Q)b(k) (58) 

R K 

The operator Q commutes with the total angular momentum operator, and the term 

A A 

im(v k • h x V) is simply the component of the angular momentum in the direction of h. 
Thus, the operator T commutes with the total angular momentum, implying that the 
functions b are indexed by l. Further consideration shows that the magnetic eigen- 
functions are simply the spherical harmonics Y^ 1 . Transforming to a new set of 
spherical coordinates (kp, 0, $), with the new polar axis along the direction of h, re- 
sults in the following simple form for equation (58): 

db 7 

i — - = r (q - Q)b. 
d$ L 
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or, 


i A Y m (0, *) = r(q - Q)Yj n (0, 4>) (59) 

L L 


The new functions yP( 0, <£) are linear combinations of the yP 1 (9, cp) with the same 
m t 

Z -value. Thus Y^ (0, $) is an eigenfunction of Q and 


m 

T Z 


-m 

q(i - oij) 


(60) 


The relation between the scattering eigenvalue problem and the magnetic scattering 
problem is clear. The presence of the vector h in the operator T removes the azi- 
muthal degeneracy present in the scattering eigenfunctions. 

Now that the magnetic eigenfunctions have been found analytically, we could proceed 
by using equation (26) for the conductivity. Although it is not necessary, it is illumi- 
nating to proceed with the calculation of the C matrix. Because of the simple form of 

■A. A __ 

T in an h-fixed set of polar coordinates, it is easier to compute the matrix h • C di- 
rectly. From equation (33), 

h • C = - ^ [(1 - aj)(l - a^ 2 / Y z m ( S )£ • X V)Y z r ( S)dS (61) 


where the symbol S has been used for (9, cp). Again, in terms of the variables (0, <3>), 
the quantity (m/h)(h • v, x V) is simply d/d$. The transformation from the spherical 
harmonics Y™(S) to the new coordinate system, resulting in Y* (0,$), is obtained 
from the fact that the yP are basis functions for the irreducible representations of the 
rotation group. Thus, from reference 34, 


Y”(S) =£ Yf' (e, *)D® n ,(h) 


m 1 


a relation that is considered in more detail later in this section. Finally, we have 


fa - ~ t 1 ~ a i> Yj MD Mm ^ D Mm ' 6 ZZ ' 


(62) 


M 
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The major difference between the present example and the discussion of the section 
SYMMETRY PROPERTIES OF THE SCATTERING EIGENFUNCTIONS arises from the 
fact that different irreducible representations of the rotation group are not coupled in 
h • C ; the matrix is block diagonal in l. The symmetry breaking provided by the re- 
striction to be a point group makes the structure of h • C more complicated. 

Summarizing the magnetic scattering eigenvalue problem, the magnetic eigenvalues 

are 


m 

T l 


-m 

q(l - op 


m = -l, -l + 1, . . . ,1 - 1,1 


and the corresponding eigenfunctions are 


hm® = Y ; m < e > *> = T. Y 7 m '®D® m .W 


i 


mm T 


(63) 


m T 


A comparison of equation (63) with equation (30) shows the correspondence between 


D mm’ and nn- 


The conductivity tensor is obtained from equation (26), with a factor of 2 included to 
count both spin states, 


a 


v \ u 



b Zm 6 ( e k ' 


e„)dk 


l, m 


(64) 


where 


(q - Q)"\ 


I C b *m 


( 65 ) 


to make the correspondence with equation (25). In the case of a spherical surface the 
velocity vector is everywhere radial, and 


mv ju = Kk p = Rk F k p 


33 



or, 


v = 



defining the unit radius vector k, with components k^ . The k^ have the property that 
they are linear combinations of only l - 1 spherical harmonics. Thus, in both equa- 
tions (64) and (65), the summations over l are restricted to the terms for 1 = 1. The 
final form for the conductivity is 


V[l 



( 66 ) 


The factor in equation (66) that plays a role similar to the g^ n of equation (36) is 
the integral 


I 

vm 


'/Vim 


(0, <f>)sin 8 d 8 d cp 


with bj m given in terms of the Y™(S) by equation (63). In order to evaluate these in- 
tegrals, a specific representation for the D- matrix must be chosen. The D is the 
3x3 (for 1 = 1) rotation matrix that describes the rotation from the fixed coordinate 
system to the h-fixed coordinate system. By convention, this is a rotation matrix for 
the Euler angles (<f>,0,-<$), and the form of D is (ref. 34) 


- (1 + cos 0) 
2 


— — sin 0e i( ^ - (1 - cos 0)e 2i< ^ 

2 


dW(*,0,-$) 


— — sin 0e -i ^ 

f* 


cos 0 


— — sin ©e^ 


— (1 - cos 0)e -2i<& 
2 


- — — sinOe - *^ 

V* 


— (1 + cos 0) 
2 
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1 m 

With this form for D, the /3* and the integrals I ym may be evaluated. The results 
are 



and 


Vm = [ D -V, m* + D !!L«] 



The magnetic scattering eigenvalues for l = 1 are obtained from equation (60) and, with 
the definition t~^ = q(l - ckj), are 


±1 

T 1 


= =FT 


0 

T 1 


= 0 


The conductivity may be evaluated directly, but the form of a itself is not of much in- 
terest. For example, from equation (66) the zz-component of a is 
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where 


zz 


e 2 r 


m 


2V 4tt 
(2tt) 3 3 


D (l) 

U 00 



1 - it or 1 + it or 


ne 2 r 

m 


2 

cos © 


2 

sin Q 
1 + (wr) 2 


n = 


2V 4n k 
(2tt) 3 3 


3 

F 


is the density of electrons in the metal. 

The more interesting result, the quantity that would be measured experimentally, is 
the resistivity tensor p^ = (cr ) . With the definitions of y = <x>t and = m/ne r, 

the resistivity tensor is 


1 


-y cos 0 y sin © sin $ 


£M = 

p 0 


y cos 0 


1 


-y sin © cos $ 


|_-y sin 0 sin $ y sin 0 cos $ 1 j 

This result is the familiar free-electron result, with the difference that the coordinate 
system has been rotated to axes in which the magnetic field is in the direction (©,$). 

If the model of a spherical Fermi surface is to have any validity for the alkali 
metals, in particular potassium, it is clear that a spherically symmetric scattering 
model is not sufficient to explain the experimental results (ref. 44). One possible model 
for the scattering, the degenerate kernel model, is the subject of the next section. 


DEGENERATE KERNEL THEORY 

The Boltzmann equation (BE) provides the framework for the description of those 
transport properties of metals that do not arise from strictly quantum- mechanical ef- 
fects such as the Landau level quantization. The nature of the BE, an integrodifferential 
equation for the distribution function, has always hampered attempts to provide calcula- 
tions for any but the simplest models for metals. The difficult part of any calculation 
has always been the manner in which the scattering term of the BE has been treated. 
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The relaxation-time approximation (the assumption of either a constant or a scalar func- 
tion of wave number to replace the effect of the scattering) has been the most common 
approximation. This achieves the desired reduction of the BE, but at the cost of aver- 
aging the scattering over all states. Thus, the knowledge about the relation between the 
topology of the Fermi surface and the galvanomagnetic properties of the metal has been 
gained at the expense of an understanding of the electron's interaction with its environ- 
ment in the metal. 

Sondheimer (ref. 21) formulated a way to treat the scattering in a more complete 
manner by making the ad hoc assumption that the kernel, Q(k, k') of equation (18), is de- 
generate. Written explicitly, this assumption is that 

N 

Q(k,k>) = b..q.(k)q.(k') (67) 

h 

j=l 

where {q^j is a set of N arbitrary functions of integrable square (ref. 24). The impetus 
for this assumption is the existence of well-documented methods of treating Fredholm 
integral equations having degenerate kernels (ref. 22). Since the summation in equa- 
tion (67) is over a finite number of terms, the integral equation reduces to a set of linear 
equations. 

The approximation of the scattering kernel by a degenerate kernel may be partly 
justified as a scattering model. Further, the basis functions q^ are not arbitrary. 

Mann (refs. 2 5 to 27) has demonstrated that, for an electron gas that interacts with a 
static system of impurities, the transition probability for scattering could be written as 
in equation (67), where the {q^j are Wannier functions for the electrons. Further work 
by Mann (ref. 27) and Seeger et al. (ref. 45) discusses the transport properties of such 
a system as a model for a metal, but as yet no treatment for the galvanomagnetic prop- 
erties of the model has appeared. 

In this section the galvanomagentic properties that result from a spherical Fermi 
surface and a simple degenerate kernel are discussed. This discussion serves three 
functions. First, an inconsistency in the Jones and Sondheimer calculation (ref. 24) is 
examined and its origin in the arbitrary truncation of the degenerate kernel is demon- 
strated. Second, the natural representation provided by the symmetry-adapted functions 
becomes even more apparent when applied to a degenerate kernel. Last, the sample 
calculations show the strong effect upon the galvanomagnetic properties for a simple an- 
isotropic scattering model, for even a spherical Fermi surface. 

It should be apparent from the discussion of scattering eigenfunctions, and in parti- 
cular the appearance of equation (45), that the degenerate kernel theory has even wider 
applicability than has been realized previously. With the assumption of the symmetry 


of the kernel with respect to the point group operations, that is. 


Q(Rk, Rk') = Q(k,'k') 

there exists a biharmonic expansion of the function Q in the symmetry-adapted func- 
tions. If this expansion is truncated after N terms, then the result is equivalent to 
equation (67). The advantage of this formulation is that the integrals of the symmetry- 
adapted functions may easily be calculated from the representation theory of the point 
group, as has been demonstrated here. The main disadvantage of working with a degen- 
erate kernel lies in the uncertain relation between the postulated kernel and the actual 
scattering processes. To make this connection satisfactorily requires the calculation 
of Q(k, k T ) from an interaction Hamiltonian, and then the use of the orthogonality prop- 
erties of the symmetry-adapted functions to invert equation (45) for the expansion coeffi- 
cients A^,. This inversion adds a substantial amount of labor to the computations - 
labor that is unnecessary in view of the present formalism which only requires the func- 
tion Q. 

To return again to the work of Jones and Sondheimer, the use of the (qj as basis 
functions allows a solution to the linearized BE for the conductivity and, hence, the re- 
sistivity. Since the basis functions are unspecified, the resulting expressions are for- 
mal in nature, equivalent to equation (23). To make possible a comparison between the 
results of reference 24 and the present work, consider the model consisting of a spheri- 
cal Fermi surface and the scattering kernel 


Q(k, k') = a Q + £ kk^ ♦ a 4 (k 4 + k; 4 ) 
a L 


+ b,k 2 k' 2 + c.(k k' 3 + k 3 k' ] 
4 a a 41 a « a a 


+ X e 4 k ot k (3 k a k h ( 68 ) 

a f^/3 

as given by equation (3. 4) of reference 24. The constants aQ, a^, b^, Cg, c^, and e^ 
must be chosen consistent with the requirements that Q(k, k T ) be normalized and 
positive-definite for all k and k T . In terms of these partially adjustable parameters, 
Jones and Sondheimer show in their figures 1 to 3 (ref. 24) the variation of the longitu- 
dinal magnetoresistance as the magnetic field varies from (100) to (111) . More im- 
portant than the quantitative results is the observation that the behavior is not that of a 
free-electron gas. There is a nonzero longitudinal magnetoresistance that initially in- 
creases but saturates as the magnetic field becomes large. Thus, a reasonably simple 
anisotropy in the scattering kernel produces a marked effect in the galvanomagnetic 
properties, even for spherical Fermi surfaces. 
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Rather than simply taking the kernel of equation (68) with some coefficients chosen 
for comparison with those of reference 24, consider in a little more detail the problems 
involved in working with a degenerate kernel. This is done by using equation (68) as a 
working example. It is easily seen that most of equation (68) is extraneous. Jones and 
Sondheimer show that only that part of Q which is odd under coordinate inversion can 
contribute to the transport current. They use this property to eliminate the problem 
of the singularity that is in the operator q - Q. This report shows that a condition 
stronger than the parity requirement may be formulated. Of all the irreducible repre- 
sentations contained in Q, only those which are not excluded by the Kronecker product 
selection rules for the formation of the C matrix need be considered. Since, for the 
cubic group, this excludes all irreducible representations with even parity, the condi- 
tion of Jones and Sondheimer (ref. 24) is included. 

All these observations may be illustrated by considering as an example the kernel 


Q(k,k’) 




r _ r 
(k)X r 


15’ 


r _ 
(k’) 


(69) 


The constant term serves only to normalize J'J' Q(k, k')dk dk' and is not important 
otherwise. The function q(k), obtained by integrating equation (69) over dk’ is simply 
1/477. This kernel differs from equation (68) in two respects, one minor and one quite 
substantial. The quartic factors multiplying a 4 and b 4 of equation (68) give integrals 
that do not vanish when q(k) is formed. This minor inconvenience is due to the use of 
nonorthogonal basis functions. This restriction of attention to a constant q(k) is of only 
minor consequence. The major difference can be seen by a consideration of the criteria 
for the choice of the coefficients, in both equations (68) and (69). The condition that 
must be satisfied, once the form of Q is chosen to satisfy the group symmetry require- 
ments, is that Q must be positive-definite. In the Jones and Sondheimer discussion 
(ref. 24) of the degenerate kernel, this meant the evaluation of Q for sufficient values 
of k and k' to give inequalities that must be satisfied by the coefficients. Equation (49) 
implies that the scattering eigenvalues are obtained from a secular determinant of the 
A^ t of equation (69). (The difference between Q and the reduced kernel 


Q(k, k')/ |/q(k)q(k') is a constant multiplier, so that eq. (69) also represents the reduced 
kernel. ) 

Explicity, since the matrix A is Hermitian (real symmetric when the real repre- 
sentaions for the Kubic harmonics are used), the secular equation is 
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A 11 - a 

A 13 

A 13 

A 33 


( 70 ) 


At this point, the major distinction between the two kernels (eqs. (68) and (69)) is ap- 
parent. The coefficients of A gg, arising from 1 = 3 spherical harmonics, are of the 
sixth degree in k and k\ The Jones and Sondheimer kernel was truncated to include 
terms only to the fourth degree. This omission of the Agg term from equation (68) has 
interesting consequences, since it implies that the scattering operator cannot be positive- 
definite. This observation may be made from equation (70) by setting Agg - 0 and, 
thus, for any nonzero A^j and A^g there is one root that must be negative. Since a 
necessary and sufficient condition for an operator to be positive-definite includes the re- 
quirement that all eigenvalues be positive (ref. 23), the truncation of the series must 
be performed carefully. 

The galvanomagnetic properties of a system with a spherical Fermi surface and a 
scattering kernel of the form given by equation (69) are easily calculated with the present 
formalism. This discussion is not meant to provide a comparison with any experiment, 
but, rather, is an investigation of the properties of a simple scattering kernel. The 
method that was chosen to accomplish this investigation was to choose a value for A^ 
and Agg (in this case the values chosen were 0. 5 and 0. 25, respectively) and to deter- 
mine the range of Ajg that gave positive, nonzero eigenvalues a. The solution of 
equation (70) for the eigenvalues and the eigenvectors, providing the appropriate coeffi- 
cients for the linear combinations of the Xj and Xg to form the scattering eigenfunc- 
tions, then provides sufficient information to compute the resistivity tensor as a func- 
tion of the magnetic field strength and orientation. 

A convenient way to present the results of the computations is to define a quantity 
that parametrizes the strength of the mixing term A^g. One such parameter is the 
ratio of the off-diagonal part to the Euclidean norm of the matrix 



The coefficients are then normalized to give a unit norm, so that the quantity A^g is 
itself the desired anisotropy parameter. The components of the resistivity tensor that 
are of interest are the longitudinal resistance p zz , the Hall resistance p xy , and the 
average transverse resistance p, = (p + p )/ 2. These quantities are calculated for 

tl AA V V a 

different values of A^g, as functions of H, and for different field orientations h. 
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The values chosen are intended to be representative of both a weak- scattering case 
(Ajg = 0.06275) and a strong- scattering case (A^ = 0.26726). This characterization of 
"weak” and "strong” is justified by the extreme dependence of the saturation magneto- 
resistance values upon the anisotropy parameter. The constant q“* is used as an ef- 
fective relaxation time r or simply r, in agreement with the usual experimental 
convention. Figures 1(a) and 2(a) show the change in longitudinal and transverse re- 
sistance as a function of c or, for fields in the two symmetry directions, (100) and 
(111) . Figures 1(b) and 2(b) give the Hall resistance p xy and the Hall constant 
(dp xy /d(o)T)), for the applied field along the (111) direction. The pairs of curves for 
each value of Ajg are qualitatively the same, but the scale has been changed for the 
larger Ajg, demonstrating the strong dependence upon the scattering. The magnetore- 
sistance curves show the initial rise (probably quadratic) and ultimate saturation that 
are expected for a spherical Fermi surface model (ref. 3). The Hall properties reflect 
this low-field behavior. Although the Hall resistance p xy is predominantly linear, 
there is a small region near H = 0 that shows a nonlinear behavior. This is indicated 
by the variation of the Hall constant in figures 1(b) and 2(b). The primary difference 
between the two scattering cases lies in the scale of the effects. As the off-diagonal 
part, Ajg, increases, the nonlinear behavior near the origin becomes more pronounced, 
saturation of the magnetoresistance requires a larger field, and the value of the satura- 
tion magnetoresistance increases quite markedly. 

To complete the discussion of this scattering model, figures 3(a) and (b) show the 
variation in the saturation values of p zz and p^ r as a function of field orientation 
around the triangle ((100) , (111), (HO), (100)). 

Although the simplified, degenerate kernel model discussed is much too simplified 
to be a serious model for a metal, the discussion of the present section shows some 
points that K ear emphasizing. There is a great deal of information contained in the 
field- independent quantities C and g . The evaluation of all components of the re- 
sistivity tensor as a function of field strength and orientation is a simple task, and the 
amount of information that can be obtained about the galvanomagnetic properties of the 
system should be apparent at this point. This complete description of the p tensor 
must be compared with previous calculations (ref. 24) that must, perforce, concentrate 
upon only one component. 

In summary, the degenerate kernal model clearly demonstrates the effect of a 
simple scattering model upon the galvanomagnetic properties of an electronic system 
with a spherical Fermi surface. Even the simple model discussed here gives results 
markedly different from free-electron behavior. This simplified model is not suggested 
as a realistic model for potassium, the metal it might be thought most nearly to repre- 
sent. On the other hand, the discussion of even this simple model demonstrates the 
power and usefulness of a method capable of treating the electron scattering. 


41 


■■■niniiiiiimiiiH ii iiiih^miiihi 


ii nn min 


CONCLUSIONS 


The independent-particle model for conduction electrons in metals and the Boltzmann 
equation have provided a great deal of information about conduction processes in metals. 
Such diverse effects as magnetoresistance, cyclotron resonance, and magneto-acoustic 
attenuation (among many others) have been explained, at least qualitatively, with the 
semiclassical theory (refs. 3 and 16). Even the descriptions of quantum- mechanical 
effects, such as the de Haas - van Alphen effect, have used the concepts provided by the 
semiclassical theory as a basis (ref. 2). The current state of experiments on metals 
indicates that the qualitative description provided by the relaxation-time approximation 
is no longer sufficient. 

If galvanomagnetic experiments are to be used to probe the interactions of the con- 
duction electrons with the environment in the metal, then clearly a theory adequate to 
support these experiments must be more detailed than any theory is at the present. The 
ideal theory would allow the inverting of galvanomagnetic data to obtain the scattering 
probability Q(k, k f ). This goal cannot be completely attained, but this report provides a 
method to partly fill the needs. The scattering probability is used in a procedure that 
gives a closed-form solution to the Boltzmann equation for a metal in static, homoge- 
neous, electric and magnetic fields. 

From a model for the Fermi surface and the scattering processes, the resistivity 
tensor may be calculated for any magnetic field strength or orientation with respect to 
the crystal axes. A feature of this formalism is the ability to separate geometrical ef- 
fects caused by the Fermi surface topology from anisotropic scattering effects. For a 
large number of metals the Fermi surface has been very accurately determined from 
band structure calculations, supported by experiment. Using this information as input 
to the procedure makes possible quantitative judgments about the effect of the scatter- 
ing. The example of the degenerate kernal model discussed in the preceding section 
demonstrates the strong role that the scattering plays in the galvanomagnetic properties 
even for a metal with the simplest Fermi surface, a sphere. The vector mean free 
path A(k) has been introduced (refs. 8 and 9) to treat as completely as possible the ani- 
sotropies inherent in the system. 

The formalism discussed in this report has been limited to the solution of the 
Boltzmann equation describing a metal in static electric and magnetic fields, for the 
case of elastic scattering. As discussed in the text, this limitation upon the scattering 
is not stringent; the formalism may be generalized to include inelastic scattering, at a 
cost of a substantial increase in the numerical work that implements the formalism. In 
addition, the linearization procedure and the method of solution of the BE should provide 
a valuable method of attack for other transport phenomena. Following the discussion of 
Taylor (refs. 8 and 9) for the thermopower of copper, this procedure is applicable to a 
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description of thermomagnetic effects in metals, currently a subject of substantial ex- 
perimental interest. 

In summary, a formal solution to the linearized Boltzmann equation, valid for any 
strength of the magnetic field, has been presented. The numerical procedure to apply 
this formal solution to any Fermi surface to investigate the scattering processes in a 
metal has been described. As an example of the power of this technique, some con- 
sequences of a scattering model introduced by Sondheimer (ref. 21) for a spherical 
Fermi surface model have been investigated. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 8, 1969, 

129-02. 
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APPENDIX A 


SCATTERING AND MAGNETIC SCATTERING EIGENVALUE PROBLEM 

In the course of solving the linearized Boltzmann equation it is necessary to con- 
sider the eigenvalues and eigenfunctions of two non-Hermitian operators. In addition to 
the existence of such eigenvalues it is also necessary to assume that these eigenfunctions 
form a complete set for the expansion of functions within a Brillouin zone. The two 
problems are, in fact, quite similar, and it is instructive to consider the simpler case 
of the field-free scattering eigenfunctions first. 

The integral operator Q is a real, bounded, positive-definite linear operator 
(ref. 23), and as such possesses a real eigenvalue spectrum. It is not necessary to as- 
sume that the kernel Q(k, k') is degenerate, so that, in general, the number of eigen- 
values is denumerably infinite, but bounded. Unfortunately, the operator which appears 
in the solution of the Boltzmann equation, as evidenced in equation (28), is q - *Q, which 
will not be Hermitian unless q is constant. Instead of this operator, the following 
eigenvalue problem can be considered 

/ - S f Y n (kW = a n Y n (k) (Al) 

fq(k)q(k') 

Since Q(k, k’) is a probability function, the integrated kernel q is a real function of k 
and will not vanish unless Q(k, k’) is identically zero for k and all k f . The kernel of 
equation (Al) is symmetric, real, and continuous, so that the existence of the eigen- 
values a' n and the eigenfunctions Y n (k) are guaranteed. 

Now define a new set of functions a n (k) such that 

Y n (k) = a n (k) f^) (A2) 

Then equation (Al) may be written as 

f Q(k,k')a n (k’)dfi' = « n a n (k) (A3) 

which is just the eigenvalue problem for the operator q *Q. Thus, the eigenvalues of 
q -1 Q are identical to the eigenvalues of equation (Al), and the eigenfunctions are obtained 
by inverting equation (A2). 
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The orthogonality condition for the scattering eigenfunctions (eq. (29)) follows 
from the symmetry of the kernel of equation (Al). By writing equation (Al) for another 
eigenfunction (e. g. , Y r (k’)) and using the symmetry of the kernel, the result is 

(“r ' “n> f Y r (k)Y n (k)dS2 = 0 

Thus, for a n * a v , the desired orthogonality is attained. If, on the other hand, 
a n = a r , the degenerate eigenfunctions may be made orthogonal by the Gram-Schmidt 
orthogonal ization process (ref. 23). From this relation then, equation (33) follows im- 
mediately by the substitution of equation (A2), that is, 

J a r (k)q(k)a n (k)df2 = 5 rn (33) 

The completeness of the scattering eigenfunctions a n can also be seen from equa- 
tion (A2). The Y n (k), being an orthogonal set of eigenfunctions of a Hermitian operator, 
form a complete set of functions. This means that, given a function f(k), there exists a 
set of coefficients C n such that 


f(k) =Yj C n Y n (k) (A4) 

n 

Since, from equation (A2), each a n is obtained from the corresponding Y by division 
by the function q(k) for all n, then 


m . 

q(k) 


E C n a nW 


(A 5) 


in which the coefficients C n are the same as in equation (A 4). Equivalently, the ex- 
pansion 


^)=E D n a n( k ) 

n 

may be inverted to obtain the D n by multiplying by the product a r q and using equa- 
tion (A3), that is, 
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(A6) 


D r = J g(k)q(k)a r (k)dJ2 

which is another way of describing completeness. 

This completes the discussion of the existence and completeness of the field- 
independent scattering eigenfunctions. Not only is this discussion important in its own 
right; it is useful in guiding the treatment of the magnetic scattering eigenvalue problem. 
The problem involved here lies in the fact that the magnetic scattering operator T, al- 
though it is the product of two Hermitian operators, is not itself Hermitian, in general. 
Consider a more general problem and write an operator T as the product of two oper- 
ators 


T = Q -1 P (A7) 

where P and Q (and thus Q - *) are each Hermitian. The operator T will not be 
Hermitian unless P and Q - * commute. There is the additional assumption that Q is 
a positive-definite, bounded, linear operator. With this assumption the formal definition 
of the square-root operator is possible (ref. 46). 

Following the discussion of the scattering operator, define a new operator T* as 
the product 


T* = Q -1 / 2 P Q“ 1//2 (A8) 

Since T T is Hermitian, the purpose of the construction, it has eigenvalues and eigen- 
functions defined by the equation 


T’x n = Q" 1/2 PQ" 1/2 x n = r n x n (A 9) 

Now, operating on the left side of equation (A9) with Q - */ 2 results in 

Q _1 p(Q" 1/2 x n ) = T n (Q" 1/2 x n ) (A10) 

Thus, a comparison of equation (A10) with (A7) shows that the eigenvalues of T are 
identical with the eigenvalues of T f (and as a consequence must also be real) and the 
eigenfunctions of T are the eigenfunctions of T' modified by the operator Q -1 ^ 2 , or 

b n (k) = Q _1 / 2 x n (k) (All) 
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Since the x n form a complete set, the same arguments as presented before indicate 
that the b n are also complete. 

This discussion applies to the magnetic scattering problem by making the associa- 
tions: 


Q - (q - Q) 

P - — (v, • h x v) 

K K 

Although this discussion is instructive when applied to the magnetic scattering operator, 
it cannot be considered rigorous. The difficulty lies in the fact that the operator P, 
being by association a derivative operator, is not bounded. There is no guarantee that 
the eigenfunctions of T lie in the space spanned by the scattering eigenfunctions a n (k). 
The question then is the existence of the expansion (eq. (34)) for the b^ in terms of the 
a n> The ultimate justification for the use of this expansion is that solutions for the ex- 
pansion coefficients y^ n can be obtained through the solution of equation (34). 
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APPENDIX B 


PROOF OF GROUP THEORETIC THEOREMS AND RESULTS 

This appendix concerns itself with the theory of group representations and the prop- 
erties of the scattering eigenfunctions of equation (28). Before starting this topic, it 
will be useful to expand upon the discussion of the basis functions of irreducible repre- 
sentations presented in the main text. 

The configuration space used here is the set of all points on the Fermi surface, 
which is considered to be a single sheet with a single-valued position vector k at each 
point. Under the action of a group S?, two points labeled k and k' are equivalent if 
k' = Rk, where R is a 3 x 3 matrix representation of a group operation St of <$ . If the 
group is of order g, then a point k will have g equivalent points (or images) on the 
surface, under Not all k will have g distinct images, and these points are usually 
referred to as ’’high- symmetry” points. 

The principal problem of interest is the effect of the group operation St upon a 
function of k. The general definition of this action is (ref. 34) 

«f(k ’) - f(k) (Bl) 

where, again, k' = Rk. The notation here is slightly deceptive. The left side of equa- 
tion (Bl) means that the value of the function Sti, whatever that is, when evaluated at the 
point k’, is equal to the value of the function f at the point k. 

The partners of a set of basis functions are defined by the equation 

*^“(k) = Yj Dj&WxM (B2) 

X=1 


■f Vi 

This equation means that, given the a? irreducible representation of <0 , there is a 
set of n ff functions, where n fl is the dimension of a, such that a transformed func- 
tion , evaluated at the point k, is a linear combination of all the basis functions. 
These functions are not uniquely determined until all the matrices of the irreducible 
representations, the D^(R), are explicitly stated. All that is guaranteed is that such 
a basis set exists for all a of &, and that there is a prescription for obtaining the 
functions, given the D matrices. 

The inverse of the group operation St must be in g?, so that equation (Bl) may be 
inverted to obtain 
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f(k’) = 9t~ X f(k) 

Thus, the basis function, evaluated at an image point is 


(B3) 

X 

since the irreducible representation may always be chosen to be unitary. The relation, 
of fundamental importance in the material to follow, relates the basis function at an 
image point to the set of basis functions evaluated at the generating point. 

It was discussed briefly in the main text that it is possible to form linear combina- 
tions of the spherical harmonics for a given l to form basis functions for some irre- 
ducible representation of a point group <& . The point of importance is that an infinite 
set of such linear combinations, in other words for all l, forms an orthonormal complete 
set of functions, spanning the same space of functions as do the spherical harmonics. 

The full discussion of the relation between the symmetry-adapted functions and the 
spherical harmonics is well documented in the literature (e. g. , the excellent article 
by Altmann and Cracknell (ref. 37) and the references quoted therein) and is beyond 
the scope of the work being described here. Let it suffice to say that there exists 
such a set of functions, characterized by X^ r (k), where 

x j r (k') = Yj *;w/s(r> (B4) 

s 

The appropriate orthonormality condition is 

f xj! r ’*(k)xj r (k)ds = 6 K ,6 jjt 6 rr , (B4a) 

The notation used, as described in the main text, is that X, r (k) is a basis function that 

th th ^ 

transforms as the r column of the j irreducible representation, and is a linear 
combination of spherical harmonics of order l . This discussion is sufficient to prove 
theorem in. 

Theorem IE: A function (k) is a basis function for the group <g if and only if 

y[(k) =^^ x j r (k) (B5) 

l 


49 



Proof: If equation (B5) is assumed to be true, then operation upon yf by an element of 
gives 

=Z V* x f « 

l 

b ,E x r^ D L(^) by e( i uation ( b2 ) 

l s 

= EE b i X l r « D Lr( R ) 

s l 

or, finally, 

s 

as desired. As usual, the ’'if” part of the proof is more involved. 

Since the X^ r form a complete set of functions, the most general expansion for a 
L tli 

function belonging to the j irreducible representation of g? is 

y, r w = Z b ; rSx l s « 

l, s 

The part of theorem III that is to be proven is that the expansion coefficients b are 
independent of the column indices r and s. This expansion, when substituted into 
equation (B2), gives 


- z b i rSx i S(k > - z b rz x i p « D i»s( R ) 

l, s l, s p 


and 

Z -Z Z hN q « D ir (R > 

r, s s Z,q 
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Equation (B2) states that these two forms must be identical. Since the X, r are ortho- 

ir* ^ 

normal, a multiplication by an and an integration will reduce both forms, and the 

result is 


E^ rSD l.s( R )=Z b ^ r (R) 

s s 

for each j, l , r, and p. The pertinent result from group theory is an orthogonality 
relation for the matrices of irreducible representations 


E^rp< R > D !,p< R >= 8 r. < B6 > 

P 

which is simply the statement of unitarity of the irreducible representations. Using this 
result then, 


b [ r =Z b rZ D ir< R > Dl pr< R ) < B ’> 

S P 

The final result is obtained from a second orthogonality relation 


Z D'' (R)D' (r) = JL 6 6 ( 

rs v ' pq v ' rp sq v 

R j 

-fVi 

where the dimension of the j representation is again nj. Thus, summing over all 
elements of the group in equation (B7) and applying equation (B8) give 


gb 


rr 

l 


-fl 

n JS 



for all r. This means that the b^ are independent of r, and completes the proof of 
theorem m. 

Theorem II (p. 23) is a statement about the form of the expansion of a kernel func- 
tion in terms of the symmetry-adapted functions and, by extension, about its eigenvalues 
and eigenfunctions. 
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Theorem II: Given a function Q(k, k’) that is invariant under all operations & of a 
group 9 (i. e. , Q(k,k’) = Q(Rk, Rk') for all dt in &), then there exists the expansion 

Q(k,ko = y 4 . X) x ] r w x ?^ *£’) 

1,1', T 

j 

Proof : The function Q(k, k’) considered successively as a function of k and a function 
of k’ certainly possesses the general expansion 

Q(k,r>=X 2 A*p ' r 'x^ r <k)x] : r ' *(k') 

i,i, v,v, 

r r* 

But Q(Rk, Rk') must have the form, from equations (Bl) and (B2), 

Q(Rk, Rk') = Y Z 4 r * j ' r ' E Xp(k)xi;P’*(k')Dj,*(R)Dj > ', pt (R) 

I,}, l',V P,P’ 

r r’ 

The last two forms must be equal; and using the orthonormality of the symmetry- 
adapted functions they result in 


Jrj'r 

hr 


'-E 


p,p’ 


A;r ! ' P ' D, rp( R ) D r'p'(R) 


A summation over the group elements and the second orthogonality relation (eq. (B8)) 
yields 


Jrj’r’ _ 1_ y A jPj ,p, 6 6 6 

IV " L U' °jj’ rr f0 pp» 


P, P’ 


The 6 -functions may be used to remove the extraneous indices, and then 


A jr . 1 PiJP 

H' ~~ L u' 

) P 
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for all r. As in the case of the proof of theorem in, it is clear that A^, cannot depend 
upon the column index r, and then the final result is the statement 

Q(k,k’) = ^ A^Xr^V) (B9) 

1 , 1 ', r 

j 


that was to be proven. 

The expansion for the kernel Q(k, k') may be integrated to obtain the function q(k). 
For the case of a spherical Fermi surface, in which the factor of the magnitude of the 
Fermi velocity in d £2 is a constant and may be removed from the integral, the orthogo- 
nality of the symmetry-adapted functions leaves only the l - 0 term for V , that is, 


q(k) 


=X A ' r ° lx * ri(k) 


(BIO) 


For a nonspherical surface, the magnitude of the Fermi velocity is not a constant, 
although it is a scalar. Thus, the integral over k’ involves an integral of over a 
linear combination of scalar symmetry-adapted functions, and the result is clearly still 
a scalar in this case. 
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Figure 3. - Saturation magnetoresistance for weak- and 
strong-scattering cases. 
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